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ABSTRACT 
High resolution X-ray imaging offers a unique opportunity to probe the nature of dust in the z < 2 
universe. Dust grains 0.1 — 1 /xm in size will scatter soft X-rays, producing a diffuse "halo" image 
around an X-ray point source, with a brightness '^ few % confined to an arcminute-sized region. We 
derive the formulae for scattering in a cosmological context and calculate the surface brightness of the 
scattering halo due to (i) an IGM uniformly enriched (fid ~ 10~^) by a power-law distribution of grain 
sizes, and (ii) a DLA-type (TVh ^ 10^^ cm~^) dust screen at cosmological distances. The morphology 
of the surface brightness profile can distinguish between the two scenarios above, place size constraints 
on dusty clumps, and constrain the homogeneity of the IGM. Thus X-ray scattering can gauge the 
relative contribution of the first stars, dwarf galaxies, and galactic outflows to the cosmic metallicity 
budget and cosmic history of dust. We show that, because the amount of intergalactic scattering is 
overestimated for photon energies < 1 keV, the non-detection of an X-ray scattering halo by Petric et 
al. (2006) is consistent with 'grey' intergalactic dust grains {Q,d ^ 10^^) when the data is restricted 
to the 1-8 keV band. We also calculate the systematic offset in magnitude, 5m ~ 0.01, for such a 
population of graphite grains, which would affect the type of supernova survey ideal for measuring 
dark energy parameters within ~ 1% precision. 
Subject headings: ISM: dust, Intergalactic medium. Large-scale structure of Universe, Scattering 



1. INTRODUCTION 

The existence of a population of intergalactic dust 
grains would have wide-reaching implications on our un- 
derstanding of the universe. In particular a "grey" popu- 
lation of dust grains, which extinct uniformly across opti- 
cal wavelengths, would affect extragalactic surveys that 
require precise optical and infrared photometry. X-ray 
astronomy provides a unique opportunity to detect and 
characterize dust that may be missed by traditional de- 
tection methods, which consider background object col- 
ors fe.g. IWrightllTMll : iMgimrd et al.ll2010B . The scat- 
tering cross-section of dust to X-ray light is highly sen- 
sitive to the grain size (of radius a), with cTgca oc a'^E~'^, 
and occurs over small angles < 1' (a//Lim)~^ {E/keV)~^ 
(|Mauche fc Gorensteinl[T98l ISmith fc Dweklll998l) . The 
scattering effect pro duces a diffuse "halo" iinage ar ound 
a point-source (e.g. I0verbec9 119651 : iMartinI [19701 ). po- 
tentially capable of being observed with high-resolution 
optics such as those on Chandra. X-ray scattering ha- 
los have been observed around bright point-sources and 
gamma-ray burst afterglows that occur behind dusty 
cloud s in the Milky Way interste llar medi um (ISM) (e.g . 
Roig[l9 83: Mauchc fc Gorenstein 1986; Witt et al.|[200ll: 
Smith et al.i i2006i: iVaughan et'al.i i2006l: iTiengo et al.l 
2010f ). Bv looking at extragalactic X-ray point sources 
(quasars) well above the Galactic plane, observers can 
characterize the quantity of dust in the intergalactic 
medium (IGM) ([Flvans et al.lll985D . We wiU derive the 
formulae for scattering in a cosmological context and dis- 
cuss interesting topics for which X-ray scattering can 
yield insight into the high-z universe. 

Stars create elements heavier than helium and thereby 
are the main contributors of dust building-blocks. The 
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mass of intergalactic dust thus depends on the cosmic 
history of star formation and dust dispersal mechanisms. 
lAguirrd (|1999n estimated that about 50% of metals at 
z ~ 0.5 are present in the IGM, and if ^ 50% of those 
metals are locked up in dust, then ri^^^ ^ 10"^. More 
recent determination of cosmological parameters would 
only change these estimates ~ 25%. Constraints on the 
amount of cosmic dust using the thermal history of the 
IGM agreed t hat »^^^(z = 3) < l O "'^ for grain sizes 
g > .01 Atm (jlnoue fc Kamavall2003D . iDiikstra fc Loebl 
(j2009f) put an upper limit on fl]^^ from 7 x 10"^ to 
1.5 X 10~*, under the hypothesis that 10% of the soft X- 
ray background is the result of dust scattered light from 
AGN. From this assumption they also showed that the 
optical depth to X-rays must be<0.15atz = 1.5. A 
previous search for a cosmolo gical scattering halo around 
QSO 1508-1-5714 at z = 4.3 bv lPetric et all l|200l yielded 
a null result. They used this to place the constraint 
^du^ < 2 X 10~^. We will show that our models are con- 
sistent with their observation and that the limits placed 
on the amount of intergalactic dust can be relaxed (i}7]). 
In Section [5] we examine the potential enrichment 
sources and the relative amounts of dust we expect to find 
in the IGM. The first stars in the early universe (Pop III) 
and feedback from early star forming galaxies would have 
enriched the IGM with metals, creating a fairly uniform 
dust distribution that is comoving with cosmological ex- 
pansion. We review the dust scattering cross-section for 
X-ray light in Section [3] and use it to discuss the case 
of scattering through a constant comoving number den- 
sity of dust grains in Section [4l Quasar absorption sys- 
tems, caused by dense regions of neutral hydrogen - often 
with traces of ionized metals - already contain signa- 
tures of dust extinction (jMenard et al.l [20051 : lYork et aTl 
I2006f ). Of individual dust sources. Damped Lyman-a 



systems (DLAs) are the densest absorbers and thus have 
the highest probabihty of producing an X-ray scatter- 
ing signature. We evaluate this possibihty in Section [5] 
by calculating the X-ray scattering profile from an in- 
finitely large screen of dust with column densities typical 
of DLAs. 

Due to the large distances associated with cosmologi- 
cal scattering, the angular extent of an X-ray scattered 
halo can put limits on the uniformity of the IGM or the 
size of a dusty cloud. Due to the nature of small angle 
scattering, the halo image is sensitive to dust at redshifts 
z < 2 and is relatively insensitive to the redshift of the 
background point source. We discuss the distances and 
timescales associated with cosmological X-ray scattering 
in Section ini Observing a scattering halo, in addition to 
finding extragalactic reservoirs of dust, can provide in- 
formation about the nature of star formation and galac- 
tic structure at high -z. In Sectio n [71 we reevaluate QSO 
1508-^4714, used bv iPetric et all (|2006l) , in the 1-8 keV 
band to show that large dust grains with flj^^ — 10~^ 
cannot be ruled out. We summarize the variety of impli- 
cations in Section [SI 

2. POTENTIAL SOURCES OF INTERGALACTIC DUST 

Intergalactic dust may be distributed uniformly as a 
result of star formation in the early universe. If that 
is the case, an azimuthally symmetric X-ray scattering 
halo may be observed around a background point-source 
at a distance large enough to provide sufhcient column 
density for scattering. Dust may also be efficiently dis- 
tributed into galactic halos or the IGM through pressure- 
driven feedback from galaxies. This effect may be par- 
ticularly important at z ~ 2 — 3, during the epoch of 
star formation, and for small galactic halos from which 
outflows can more easily escape. Previous work has con- 
cluded that dust grains > 0.1 /xm are more efficiently ex- 
pelled than smaller grains because they are grey to opti- 
cal light, hence receiving more radiation pre ssure, and are 
less susceptible to decelera tion by gas drag (jDavies et al.l 
119981: iFerraraet a'Llll99H) . Finally, if very little dust 
(r^Ji^^ < 10~^) is evenly distributed throughout the cos- 
mos today, we could potentially observe X-ray scattering 
from objects that are already known to have high column 
densities of gas and metals. Quasar absorption systems 
offer a convenient opportunity to observe X-ray scatter- 
ing from dust that may be present in these gas reservoirs. 

2.1. The Earliest Star Formation 

Population III stars formed from the near-pristine pri- 
mordial hydrogen gas originating from the Big Bang, are 
typically '^ 40 — lOOM©, and fo rm in dark matter halos 
~ lO^Mtn between z '^ 20 - 30 dBromm fc Larsonll2004l: 
iHosokawa et al.ll201lD . These stars are very large be- 
cause cooling below 10^ K can only proceed by collisional 
excitation of the rotational modes of H2 molecules. Frag- 
mentation into smaller clumps is possible once gas be- 
comes enriched to a particular threshold, above which the 
cooling rate from metals dominates over H2 cooling and 
exceeds the heating rate due to gra vitational collapse: 
[C/H 1 « -3.5 and [0/H] w -3.1 (IBromm fc LarsonI 
120041) . Many consider this critical metallicity to mark 
the transition between the era of Pop III and Pop II 
stars, which have low metallicity but are able to form 



smaller (^ 1 Mq) stars that are observed today in glob- 
ular clusters, dwarf spheroidal galaxies, and galactic ha- 
los. The large characteristic mass-scale of Pop III stars 
indicates that the majority of stars ended their lives as 
supernovae, distributing metals. If some of these met- 
als are locked up in dust, it would contribute to cooling 
and allow gas to fragment at metallicitie s well below the 
critical threshold for metal line cooling ([Schneider et al.l 
l2006f ). The recent observation of a l ow-mass, extremely 
metal poor star with [Z/H] -- -5 ([Caffau et al.l 120111 1 
supports the hypothesis that the first generation of stars 
produced dust. 

Because Pop III stars are massive and form in small 
halos, feedback from radiation and thermal injection 
from supernovae are significant. Thus the first gener- 
ation of stars can efficiently rem ove gas from mini-halos 
and mix metals into th e IGM ( Hcger & Wooslcv 2002[; 
IBromm et aD[2003l: [Abel et a l. 2007:' Wise & Abel 2008). 
In simulations bv lWise fc Ab el (2008), the second gener- 
ation of Pop III stars had metallicities —5 < [Z/H] < —3 
and formed around z = 20; the third generation reached 
metallicities [Z/H] > —3. These serve as an upper limit, 
because the simulations assumed that every Pop HI star 
ended in a pair-instability supernova (PISN). At solar 
metalli city, the mass ratio of metals versus H and He is 
~ 0.02 (|Drainell2011[ ). Using this value to scale the metal 
mass by metallicity, the total dust-to-gas mass ratio of 
the Pop HI enriched IGM is 
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where /^cp is the mass fraction of metals locked in dust. 
Dust models suggest that /dep can be rat her high for both 
PISNe (0.3-0.7) and Type-II S Ne (0.2-1) (ISchneider et al.l 
12004 iTodinifc Ferraral I2001D . Taking the metalficity 
threshold above which hydrogen clouds fragment, we can 
make a reasonable guess that the first stars may have en- 
riched the IGM to a uniform background level of [Z/H]^ 
-4. Equation [U then implies (Md/AfH)iii < 2 x 10"^. 

Compare this to the ratio implied by assuming a par- 
ticular value for the dust mass density in the IGM: 



f 






7 X 10"'' n 



IGM --1 
dust "-H 



hh 



(2) 



IGM 



where /175 is Hq in units of 75 km/s/Mpc, and nn ~ 
10~^ cm""^ is the average baryon density of the universe 
today. If f^du^ ^ lO"'^, the first generation of stars 
would have contributed ~ 1/350 of that dust. If Pop 
HI stars were responsible for the entirety of IGM dust, 
they would have to pollute the universe to [Z/H]~ —1.5, 
and all of those metals must be locked up in dust. This 
is not likely to occur before the onset of smaller stellar 
populations and galaxy evolution. Feedback from the 
second generation of stars and early galaxies is therefore 
required to bring the overall dust density of the IGM to 
a threshold detectable via X-ray scattering. 

2.2. Feedback from Galaxies at High and Low-z 

There is direct evidence that L* galaxies today can ex- 
pel dust out to a few kpc at least. Polarized light, a 
consequence of scattering, was observed around the disk 
of the starburst galaxy M82 well before its dust y out- 
fiows were imaged by Hubble and Spitzer (,Schmidt et al.l 



119761) . The outflows are fueled by starburst activity and 
can host super-solar abundances, a result of supernovae 
enri chment, effectively c ontaminating the IGM with met- 
als (|Konami et alJl2011[ ). 

In theory, dust grains can be efficiently expelled from a 
galactic disk when the forces of radiation pressure, grav- 
ity, and gas drag are weighed. Large grains (> 0.1 ^m) 
can be expelled completely from optically thin disks, 
while small grains are either retained close to their for- 
mation sites or trapped in layers within the galactic halo 
(JDavics ct al. 1998; Grccnbcrg ct al. 1987). The expul- 
sion is also more efficient for graphite grains ^0.1 /im, 
which absorb more light across the visual spectrum than 
silicate grains of comparable size, and can reach veloc- 
ities of several hundred km/s. The smaller timescale 
for exiting the halo allows graphite grains > 0.1 /im to 
survive sputtering processes compared to silicate grains, 
which exit slowly and are thus more likely to be destroyed 
(jBarsella et all 119891: iFerrara et aIlll99H) . On the other 
hand, graphite grains can accumulate more charge than 
silicate grains as they exit a galactic disk. Coulomb forces 
prevent graphite grains from traveling as far as silicate 
grains with the same starting velocity, an effect that may 
allow silicate grains to popula te a galactic halo or the 
IGM ()Bianchi fc Ferraral 120051) . Regardless, enrichment 
processes suggest that an IGM dust population is likely 
to contain grains > 0.1 /im in size and is consequently 
grey to optical light. 

Dust grain expulsion is not just a side-e ffect, but may 
be nec essary for galactic outflows to occur. iMurrav et al.l 
()2005[) argued that the Faber- Jackson relation between 
an elliptical galaxy's luminosity and velocity dispersion 
can be explained by the balance between momentum 
driven winds and stellar luminosity. They also showed 
that, for radiation pressure to act efficiently, dust is re- 
quired to contribute to the total ISM opacity. As mate- 
rial is expelled, the gas may become sufficiently diffuse 
for the dust to decouple and migrate outwards alone. 
They estimated the radius at which this occurs [R^g), 
with respect to the velocity dispersion (ct), gas fraction 
(/g), and dust grain size (a): 
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where 0.1 /im, 200 km/s, and 0.1 arc fiducial values for 
a, CT, and fg respectively. Dust grains ^ 0.1 /im expelled 
from very small halo potentials (ct ~ 50 km/s) would 
decouple from the gas 10 kpc beyond the site of origin. 
Thus during the epoch of star formation, gas may not re- 
liably trace the spatial distribution of intergalactic dust. 
If that is the case, the fraction of IGM metals locked up 
in dust may be rather large {'^ 0.5). 

An estimate for the total metallicity of the universe 
as a function of redshift can be obtained by integrating 
over the observed star formations r ates, yielding Q,z = 
2.5 X 10"^ at z = 2.5 (|Pettinill200"6[) . The total mass of 
metals observed in Ly-a absorption systems and Lyman 
Break Galaxies (LBGs) aX z = 2.5 differs by the amount 
predicted by a factor of few, co ined the "missing metals 
problem" (|Pageil l200l iPettinil [2004) . The discrepancy 
may be due to uncertainty in the star formation rates 
between z = 4 — 10 or galax ie s that have be e n miss ed by 
UV-selection (jPettinil |2006| ) . iBouche et all (|2007D used 
analytic calculations to conclude that around 50% of the 



total expected metal mass at z = 2 could have been 
ejected by sub-L* (the majority of which comes from L 
< 1/3 L*) galaxies. If 50% of intergalactic metals are 
locked up in dust, then rjJi™(z < 2) - 8 x 10"^ due 
to pollution by sub-L* galaxies during the epoch of star 
formation. X-ray scattering halos can help test whether 
or not metals really have been missing from observations 
of the z = 2 universe, and whether feedback mechanisms 
can account for the supposed deficiency. 



2.3. Damped Lyman-a Absorbers 

DLAs are dense regions of neutral hydrogen (TVhi > 
^) observed via broad Lyman-a absorption 
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features in quasar spectra. They account for the ma- 
jority of neutral hydrogen out to z = 5, making them 
likely candidates as the gas reservoirs for star forma- 
tion. Mg II and other singly ionized metals are also 
found at redshifts coincident with DLAs - enrichment 
that supports the link to star formation. Past surveys 
have used Zn and other volatile elements (Si, S, and 
O), showing that the majority of DLAs have sub-solar 
metallicitics. The median metallicity is [Z/Zq] ~ —1.2 
([Pettini 2006), which increases with decr easing redshift 
from [Z/H] = -2.5 to -0.5 (| Wolfe et al.l[2005r i. DLAs 
are also likely to contain dust, as the amount of Fe - 
a refractory element that condenses at low temperature 
- is observed to be depleted re lative to Zn and Si with 
increasing metal abundance (Wo lfe et al.| [2005). The ex- 
tinction profiles for Mg II and other metal absorbing sys- 
tems also show significant eviden ce for reddening by dust 
(jMenard et al.ll2005l : lYork et al.1 [20061. 

The link between DLAs and dust is strong, but the 
exact nature of these systems is unknown. If they are 
galaxies, observing their emission is difficult because they 
lie along the line of sight to quasars. Their sub-solar 
metallicities also make it difficult to reconcile the idea 
that DLAs are the progenitors to the galaxies we see 
today. They may be dwarf galaxies, inflow, or outflow 
gas associated with high-z bulges or star forming re- 
gions (Wolfe et al. 2005). X-ray scattering can test these 
hypotheses because (i) it requires a bright X-ray point 
source, which the background quasar can provide; (ii) 
the redshift to the dusty body is known, removing some 
of the degeneracy in interpreting the scattered image; 
and (iii) the dust distribution and size of the dusty re- 
gion places limits on the observed angular size of the 
scattering halo (^5]). 



3. X-RAY SCATTERING CROSS-SECTION 

The Rayleigh-Gans (RG) approximation is very often 
applied in the case of X-rays and dust because the in- 
dex of refraction is close t o one and thus the scattering 
angle s are very small (e.g. Ivan de Hulstlll98il : lOverbecM 
119651 ). The approximation treats a scattering particle 
as a collection of inflnitesimally small Rayleigh scat- 
tering regions. The phase functions for the collection 
of Rayleigh scatterers must be integrated over the vol- 
ume of the particle , and relies primarily on its geom- 
etry (jMartinl I1970D . For spherical particles with ra- 
dius a, the angular dependence for differential cross sec- 
tion can be approximated with a Gaussian distribution 



(jMauche fc GorensteidHQSGD : 
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where m is the complex index of refraction and the width 
of the scattering angle distribution is 

~ 1.04 arcmin 

'^"~F n • (^) 

Here a^m = a/ /xia and -BkcV = £^/keV. The total scat- 
tering cross-section is 

2 
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and is independent of the Gaussian approximation. We 
follow ISmith &: Dwekl ([l998) in applying the Drude ap- 
proximation, so that 
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where re is the classical electron radius and Ue is the 
number density of electrons, which we calculate from the 
mass density. For values characteristic of dust grains, 



o-RG = 6.18 X 10 
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w here p^ = p/3 g cm~'^, a typical density for dust grains. 
ISmith fc DweU ()1998D have shown that the Rayleigh- 
Gans approximation matches the more exact Mie scat- 
tering solution when a^„n <> -^kcV- For the 0.5 keV pho- 
tons and 0.25 /xm grains, the magnitude of the scatter- 
ing cross-section can differ by 20% (for graphite grains) 
to 100% (for silicate grains). There is also a limit- 
ing angle beyond which Equation |4] no longer matches 
the exact RG or Mie solution. For 2 keV photons and 
0.1 ^m (0.25 ^m) grains, the solutions diverge around 
500" (200"). For soft X-rays < 1 keV, the Gaussian ap- 
proximation matches the shape of the Mie solution out 
to at least 1000". However, absorption effects are no 
longer negligible below 1 keV, and the RG approxima- 
tion causes an overestimate for the halo brightness by a 
factor of 2 — 10 at 6'scat — 100". This effect also increases 
dramatically with grain size. For more details, we refer 
the reader to the original paper. 

< ii^kcV I'ule of thumb holds. 



Whether or not the 
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we can be reasonably confident that Equation |4] well de- 
scribes the shape of the halo image out to ~ 100" and 
brightnesses for photon energies > 1 keV. A formalism 
that incorporates Mie scattering is required to extend the 
accuracy of the calculation to energies ~ 0.5 keV. 

4. X-RAY SCATTERING THROUGH A UNIFORM IGM 

Here we derive th e integral for scattered light similar to 
lEvans et"all ()1985[ ). but for an updated ACDM cosmol- 
ogy. The path taken by light scattered at cosmological 
distances can be traced over three events: 

1. Light of frequency fom is emitted from a point 
source at redshift Zs and reaches a patch of dust 
grains at redshift Zg. 

2. A portion of the incident light is scattered by the 
dust grain patch. 




Fig. 1. — A diagram illustrating the scattering geometry in co- 
moving coordinates. Photons that scatter through angle 6scat from 
a patch of dust grains at redshift Zg will be observed an angular 
distance a away from the center of the point source. 

3. Scattered light projected onto 6'scat will reach the 
observer at z = and be observed an angular dis- 
tance a away from the center of the X-ray point 
source. 

We use the notation Uc = {1 + z)^^ to represent the 
cosmological scale factor throughout. Everywhere else, a 
denotes grain radius. 

Figure [T] illustrates the geometry of the system using 
X to represent the comoving coordinate distance: 
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We assume a flat ACDM universe such that H{z) = 

Hoy^fl„i{l + z)^ + nA, where Ho = 75 km/s/Mpc, 
r^m = 0.3, and i}\ = 0.7. 

In general, the specific fiux from a source at cosmolog- 
ical distances is 
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where R is the proper radial distance between the ob- 
server and the source. In the reference frame of a dust 
grain at Zg, the radial distance to the light source is 
Rgs = Xss/(1 + Zg). The scale factor of the universe, 
according to the dust grain, is ac,g = (1 -I- Zg)/(1 + Zs). 
The flux at the site of the dust grain is therefore 
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We use the subscript notation v, g to represents the fre- 
quency of light at the dust grain site [vg — Vemi^ + 
Zg)/{1 + Zs)] and 1^,0 for the frequency of light at the 
site of the oljserver [vq = Vem/i^ + ^s)]- 

When light interacts with the patch of dust grains, a 
small amount will be scattered, with intensity 
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where n{zg) is the number density of dust grains at the 
scattering site and dR is the small length 
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Only light scattered onto the angle ^scat = oi/ x will be 
received by an observer at z = 0, using the parameteri- 
zation X = Xgs/Xos- Light from the scattering event will 



reach an observer at z = with intensity 
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Combining Equations [TT] through [U gives the observed 
intensity of scattered hght from the patch of dust grains 

at Zr, 
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We can also substitute n{z) = ndz) (1 -I- z)^, where 
nc{z) represents the comoving number density as a func- 
tion of redshift. Drop the subscript for integration over 
all z and substitute the flux of the source as observed at 
z = 0: 
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to get the total intensity of scattered light as a function 
of angular distance from the center of a point source. 
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The differential cross-section must be evaluated at the 
frequency encountered by the dust grains [vz — vq{\+z)]. 
Equation [T71 is missing the effects of attenuation along 
the path traveled by scattered light. For a sufficiently 
homogeneous IGM, the optical depth to extinction along 
the path of scattered light (Tgs -\-Tog) will be nearly iden- 
tical to that encountered by light traveling in a straight 
line between the source and observer (tq^)- The extra 
distance traveled by scattered light, Ax ~ 10 kpc fi i6.2p . 
is not large enough to accumulate an appreciable col- 
umn density for extinction^ Under this assumption, a 
factor of 6^"^°= should modify Equation [T71 Then the 
source flux term can be replaced with the apparent flux 
'^"o e~'^°=. Thus throughout the paper we ex- 
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press the halo intensity as a fraction of the central point 
source apparent brightness. 
We use the notation 
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to describe the normalized surface brightness as a func- 
tion of angle. The total optical depth to scattering is 
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for a point-source at z = z^. Approximately r^ of the 
source photons will scatter twice, but for the cases con- 
sidered below, the optical depth to X-ray scattering is 
always Tx < 0.1. Therefore we do not need to consider 
light that might be twice-scattered into or out of the ob- 
server's sight. 

^ A notable exception occurs if a galaxy happens to lie along the 
line of sight. 



TABLE 1 

Optical depth to X-ray scattering for 1 keV photons 



Dust Model 

Power Sizes (/.tm) 


dr^/dR 
[% Gpc-i] 


z = 1 


2 3 


4 


1.0 
(p = 3) 0.1-1.0 
(p = 4) 0.1-1.0 


1.68 
0.92 
0.66 


5.2 
2.8 
2.0 


8.1 9.9 
4.5 5.5 

3.2 3.9 


11.2 
6.2 
4.4 



4.1. Dust Grain Distribution 

We use ridust — 10~^ as the flducial value of comoving 
IGM dust mass. The number density of dust grains for 
identically sized spherical particles is then 
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where h is Hq in units of 75 km/s/Mpc. The dust we 
might expect to find in the IGM would have been ex- 
pelled as a consequence of galactic feedbac k on the ISM. 
We follow the interstellar dust models of 'Mathis et al 
(|1977l hereafter MRN) andWcingartner & Drainc (2001] 
hereafter WDOl) to choose a power law distribution of 
grain sizes. 

The largest grains may be expelled from galaxies more 
efficiently because they contribute to the opacity of sur- 
rounding gas, receive more radiation pressure, and can 
withstand gas drag and sputtering processes that slow 
down or destroy small grains (i j2.2p . For the lower end 
of the size distribution, we choose 0.1 fim; dust grains 
smaller t han this will be des t royed by T ~ 10^ — 10^ K 
halo gas (jFerrara et allUQQll IAguirrelll999( ). We choose 
an upper limit of 1 /im to balance the easy expulsion 
of large dust grains with the existence, but severe drop 
in number, of grains > 1 /im. A power law with p = A 
matches the Ry — 5.5 WDOl graphite distribution in the 
region 0.1 < a < 1.0 /im. To explore the possibility that 
the dust population resulting from outflow has a larger 
proportion of '^ 1 pm size grains, we also try a p = 3 
power law throughout the paper. In the limiting case 
that all intergalactic dust is large, we will compare re- 
sults to a dust population solely comprised of 1 /im sized 
grains. 

The RG scattering cross-section as used in this study 
does not incorporate subtleties associated with different 
grain materials (silicate or graphite) . The complex index 
of refraction m is calculated from the number density of 
electrons, which for convenience we assume is uniform 
regardless of grain material (Equation ^ . Only the size 
range and power-law slope of the grain distribution affect 
the X-ray scattering intensities presented throughout 
this work. 



4.2. The Surface Brightness of Scattered X-ray Light 

As an X-ray photon travels between the distant quasar 
and the observer, its energy will shift towards the softer 
end of the spectrum. Thus the total RG scattering cross- 

^^, where Eg is the energy 
The total optical depth of 



section is o-rg oc a Eq i^+z] 

of the light observed at z = 0. 

the universe to X-ray scattering by a uniform, comoving 

dusty IGM becomes 
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Fig. 2. — The normalized surface brightness of scattered X- 
ray Ught from a point source at z = 1 (black) and z = A (grey) 
using two dust grain distributions: dn/da oc a~^ (solid line) and 
dn/da oc a~* (dashed line). The total integrated fractions for the 
p = 3 model are shown next to their respective curves. The light 
grey line illustrates the 100 ks Chandra flux limit, normalized by 
lO"'^^ erg/s/cm^ and divided by solid angle. 
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where a is the grain radius. Table [T] shows the optical 
depth for each dust grain model and Eq = 1 keV. All 
of them are significantly below the t^{z = 1) < 0.15 
hmit set by iDiikstra fc Loebl (|2009( ). For p < 5, the 
optical depth is dominated by the largest grains in the 
distribution, and the values in Table [T] scale roughly as 
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The energy dependence of the RG differential cross- 
section cancels everywhere in Equation |4] except the 
Gaussian term. To absorb the effect of cosmological red- 
shift, we use the parameter 
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because the Gaussian has a width cr oc E^^ . The scat- 
tering halo intensity as a function of observation angle 
(Equation [T7| becomes 
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For very small optical depths. Equation [23l should equal 
Tx when integrated over solid angle. 

Figure [2] shows the integrated halo profiles for 1 keV 
X-rays at z = 1 and 4. The total integrated fraction 
is shown next to each curve for the p = 3 dust grain 
model, and in each case agrees with the total optical 
depth to within a fraction of a percent. The increased 
halo intensity that results from using the p = 3 grain size 
distribution is not dramatic, around 50% brighter than 
the p = A case. This demonstrates that the halo bright- 
ness is more sensitive to the maximum grain size than 
the power-law exponent. The light grey line in Figure [2] 
corresponds to the radial profile that would arise if the 
brightness within 1" spaced annuli summed to the 100 
ks Chandra flux limit (4 x 10""'^^ erg/s/cm^), normalized 
by a point source flux 10^^^ erg/s/cm^. This illustrates 



TABLE 2 

Half-light radius for 1.0 keV scattering halo around a 

source at varying redshift. 



Dust Model 


Half-light Radius (") 
z = 1 2 3 4 


1 /im 
p = 3 
p = 4 


22 18 16 15 
31 25 22 20 
35 29 26 24 



a rough threshold at which it becomes difficult to distin- 
guish signal from background. 

The scattering halo becomes more compact with larger 
source redshift, because the characteristic scattering 
angle decreases with increasing photon energy (Equa- 
tion [5]). For a fixed observed energy band, photons from 
a high-z source scattered at a higher energy than pho- 
tons from a low-z source. Table [2] shows the half-light 
radius (containing 50% of the total scattered light) for 
1 keV X-rays. Across all grain-size models, scattering 
halos become more compact by 30% between redshift 1 
and 4. This creates a dilemma for instruments lacking 
angular resolution ^ arcseconds. Point sources at high 
redshift may have considerably larger column density of 
dust, but the ability to resolve the resulting scattered 
halo image will become more difficult. 

5. X-RAY SCATTERING THROUGH AN IGM CLUMP 

Though the universe is homogenous on a very large 
scale, observation along a single sightline is subject 
to fluctuations in baryon density, as evidenced by the 
Lyman-a forest and other absorption systems ubiquitous 
to quasar spectra. We investigate the possibility that an 
overdensity in the IGM can be observed individually via 
X-ray scattering through a screen at intermediate red- 
shift. As before. Equation [11] describes the point source 
flux from Zs onto a screen of scattering particles at Zg. 
Equations [l^] and [T31 implemented with RG scattering 
(Equations |4] and [22)) , produce 
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We use the general properties of DLAs as a fiducial ex- 
ample because they are the densest baryon reservoirs 
observed along the line of sight to known X-ray point 
sources and are thus the most likely candidates to pro- 
duce an individual scattering halo. In addition, we gauge 
which will contribute more to the brightness and shape 
of an X-ray scattering halo - the diffuse IGM or a dense 
cloud. 

The extinction curves measured so far in quasar ab- 
sorption systems best match the extinction by SMC- 
like dust. That is, there is no prominent 2175 A fea- 
ture, but smooth reddening when comparing the spec- 
tra of quasars wi t h Mg II absorbers t o tho se without 
([York et al.l [20061: IMenard et akl [2005I . [200l . Extinc- 
tion curves for the SMC can be well-approximated by 
an MRN distribution composed purely of silicate grains 
- absent are the smallest graphitic grains (PAHs) respon- 
sible for the 2175 A bump (|Peilll992D . We therefore ex- 
tend our study to include a distribution of SMC-like dust 
{p = 3.5 and 0.005 < a < 0.25 /im), appropriate for star- 
forming dwarf galaxies. A distribution of large grains 



TABLE 3 
Potential dust distributions for the average DLA 
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Name 


Dust Model 

Power Sizes (/.tm) 
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r^ 


SMC-like 

Large 

Large 


(P = 3.5) 
{p = 4) 
(P = 3) 


0.005 - 0.25 
0.1-1.0 
0.1-1.0 


3.2 X 10" 
1.7 X 10" 
6.6 X 10- 


-10 
-13 

-14 


~ 0.5% 
~2% 
~3% 



Optical depth to X-ray scattering at 1 keV for the column den- 
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10^ 



' and 2 = 0. Note that Tx 



Zg)' 



may be more appropriate for dusty regions produced by 
radiation pressure driven outflows. As we will show be- 
low, the brightness and the width of the X-ray scatter- 
ing halo offers a diagnostic for differentiating between 
the two dust populations and their relation to quasar 
absorption systems. 

The mass ratio of dust to gas in extragalactic systems 
may vary with metallicity and size distribution in com- 
parison to that observed locally. We introduce the dust- 
to-gas number ratio, k = Nd/N-a, which absorbs the size 
distribution, metallicity, and dust-to-gas mass ratio. If 
we assume that dust abundance is roughly proportional 
to metal abundance, 
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where the dust-mass to hydrogen gas-mass ratio of ref- 
erence is that of the Milky Way, {Md/MB_)MW ~ 0.009 
()Drainelf2011[ ). A relatively constant ratio of dust mass 
to metallicity is observ ed locally in studies of the Milky 
Way, LMC, and SMC (Pel '1992). We apply this finding 
liberally as a means of estimating the magnitude of X- 
ray scattering by dust reservoirs; it may not well apply 
for dust in outflows, halo gas, or the IGM. If DLAs are 
primordial galaxies, dust detection via X-ray scat tering 
provides a check on dust evolution models (e.g. Ilnoud 
[20133) . Table [3] shows the dust-to-gas number ratio and 
corresponding optical depth for several dust size distri- 
butions, assuming Z — 1/15 Zq. Note that the total 
optical depth to X-ray scattering, 

Tx ex fc E^^ (1 + Zg)-^ 

will scale proportionally with metallicity and gas-to-dust 
mass ratio. 

For 1 kcV photons and a source at z^ —2, the scattered 
surface brightness profiles from an infinitely large screen 
with the hypothesized dust column densities are shown in 
Figure [21 For comparison, the scattering halo expected 
from dust distributed uniformly throughout the universe 
is also shown (solid black). The other black lines show 
screens at z = 1 (dashed) and z = 0.5 (dotted) using the 
Large p = 4 dust distribution. The result for a screen 
of SMC-like dust is also shown (solid grey). When the 
characteristic scattering angle is smaller (e.g. with higher 
energy), the halo profile is more centrally focused. This 
explains why the surface brightness profile for the z = 1 
(Large dust) screen, despite having a smaller Tx, appears 
brighter than the z — 0.5 (Large dust) screen for obser- 
vation angles < 10". When the characteristic scattering 
angle is larger (e.g. with smaller dust grains), the major- 
ity of the scattered light is in the wings of the halo profile. 
Thus the curve from the SMC-like screen appears signif- 
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Fig. 3. — Surface brightness profiles for X-ray photons scattered 
from a DLA-type screen of dust particles and a point source placed 
at 2s = 2. The large dust distributions are shown for screens at 
z = 1.0 (dashed black line) and z = 0.5 (dotted black line). For 
angles < 10" , both profiles are flatter than that expected from X- 
ray scattering through uniformly distributed f^Jj^^ = 10~^ (solid 
black line). The SMC-like dust distribution is shown for a screen 
at z = 1 (solid grey line). The light grey vertical lines mark the 
radial size (Da) of an object at z = 1, ranging from 0.1 to 1 Mpc 
using 100 kpc sized steps. 

icantly dimmer than one might expect in relation to the 
other curves plotted in Figure [Sj 

It is evident from Figure [31 that the expected DLA 
dust column produces a scattering halo that is several 
orders of magnitude dimmer than what is expected from 
^dust ^ 10^'°^, particularly within 10" of an X-ray point 
source. Therefore, a typical DLA along the line of 
sight to an X-ray bright quasar would not be identifi- 
able with X-ray scattering unless the total dust density 
^du^ < 10"^ - 10"^. However, a high-metaUicity DLA 
(or an extreme dust-to-gas ratio) would present a signif- 
icant dust column density for scatteri ng. The rare DLA 
with Z r^ Zq (e.g. iFvnbo et al.ll2011| ) would have t^ ^ 
0.08 — 0.45 (1 -I- Zg)~'^, and double-scattering effects will 
become important. Eleven DLAs with [Zn/H] > —0.5 
(roughly correl ating with total metallicity) are present 
in a sample by iProchaska et al] f|2007. ) . suggesting Tx > 

0.02-0.14 (l + zg)-"^. 

The shape of the halo profile provides a crucial diag- 
nostic for differentiating between a uniformly distributed 
IGM and a discrete dust source. The uniform IGM case 
is equivalent to stacking many dust-scattering screens be- 
tween the source and observer. The more distant screens 
contribute only to scattering observable at very small 
angles - close to the point source - giving the radial 
profile a cusp near the center. For a scattering screen, 
the halo surface brightness has the same shape as the 
scattering cross-section (Equation[Ml) , which we approxi- 
mated with a Gaussian function (Equation^]) . Therefore, 
the radial profile of scattered light from a screen is very 
fiat for angles below the characteristic Gaussian width. 
The turn-over point, at which the brightness profile falls 
rapidly, is a function of (i) the dust grain size distribu- 
tion and (ii) the distance between the X-ray point source 
and dust screen. Dust distributions that are weighted 
towards smaller grains form broader halos because the 
differential cross-section allows for larger characteristic 
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Fig. 4. — The physical radius associated with ziimiti using a 
fixed photon energy and grain size. The diamond points mark the 
observation angle a associated with each ^umit in intervals of 0.2', 
the values increasing from right to left. For clarity, a points are 
only plotted up to 1' on the curves below Zs = 4. The results for a 
and Dyi shown here can be scaled by cro{arcmin) Ri l/{afjm£'kov)- 

scattering angles (Equation [S|). Wider halos are also ob- 
served when the light rays approaching the dusty screen 
are nearly parallel; in comparison, a screen very close to 
the X-ray point source must scatter light through larger 
angles in order to reach the observer. 

For a DLA scattering halo to be observable, the dusty 
region must subtend a large enough angular diameter to 
be viewed at sizes comparable to the angular width of the 
scattering cross-section. The angular diameter distance, 
which describes the physical size of an object viewed at 
cosmological distances. 
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is relatively constant beyond z ^ 1, for a fixed observa- 
tion angle. The light grey vertical lines in Figure [3] mark 
the angular diameter distance for an object at z = 1, in 
100 kpc steps, from 0.1 to 1 Mpc. Most objects 2: > 1 
must be larger than 200 kpc for the brightness profile 
turn-off point to be observed , but DLA sy stems are tens 
of kpc large, not hundreds (jPettinil 120061 . Observable 
X-ray scattering in the vicinity of quasar absorption sys- 
tems is likely only if these systems are associated with a 
nearby galaxy, produce outflows or are a tracer of galac- 
tic outflow themselves. If the diffuse distribution of the 
IGM is sufficiently low. X-ray scattering offers an avenue 
for exploring the structure around DLAs. 

6. REDSHIFT SENSITIVITY 

The ability for a region of dust to scatter light towards 
an observer depends primarily on the scattering angle 
and the width the Gaussian distribution in Equation 21 
We showed with Equation [22] that one can imagine all 
observable scattering is limited to angles 6cS that fall 
within a distribution of width ctq = a{Eo). Most light 

will scatter within OeS < 2(To. Given a particular grain 
size and energy, for every observation angle a there will 
be a limiting redshift (^umit) beyond which effectively 
no light will be scattered towards the observer. Smaller 
observation angles are able to probe regions closer to the 



X-ray point source, while larger observation angles are 
sensitive only to regions close to the observer. 

For 1 keV photons and 1.0 iim. grains, the effective scat- 
tering angle is limited by 2ao ~ 2'. Figure H] shows the 
angular diameter distance Da for objects as a function 
of ^limit (found by solving 0cs — 2uq numerically). The 
curves border the volume of the IGM that contributes to 
a scattering halo. Diamond points on each curve mark 
the observation angle a in 0.2' sized intervals. To extend 
the results of Figure |4] to other grains sizes and photon 
energies, the a and Da values should be scaled by ctq in 
units of arcmin. 

The contents of Figure |4] give broad insight for the abil- 
ity of X-ray scattering to probe different regions of the 
IGM. The contribution to the scattering halo for angles 
a > 0.2' and a point source at Zg = 4 covers a depth 
z < 2 of a uniformly dusty IGM. An observer who wants 
to learn more about z > 2 dust must be able to mea- 
sure the halo brightness profile for angles within ~ 10" 
of the point source center. The maximum for each curve 

in Figure H] occurs when a w ctq (1'). Therefore, the 
largest detectable extragalactic structures are at z ~ 0.5 
and have a physical diameter ~ 500 — 700 kpc, depending 
on the redshift of the background point source. 

When measuring the width of a scattering halo, there 
is ambiguity between a large object far away and a small 
object nearby. From the Zs = 4 curve in Figure IH a 
dusty cloud 200 kpc in radius would create a scattering 
halo of radius 0.2' if placed at Zg = 2. From Eauation[26l 
a 90 kpc cloud ed Zg = 1 happens to have an angular ra- 
dius of 0.2' and would thereby create a scattering halo 
of comparable size. Knowledge of the object's redshift is 
required to break the degeneracy between size and dis- 
tance, making quasars with known absorption systems 
more favorable candidates for X-ray scattering studies. 
Due to the nature of small angle scattering, halo light 
appearing at a « 2aQ will always have come from dust 
close to observer. This fact may help mitigate the effects 
of foreground dust, which will contribute to the scatter- 
ing halo several arcminutcs from the point source. 

Each annulus of an X-ray scattering halo image is sen- 
sitive to the interplay of several parameters (a, E, and 
z). We examine the curve 



d 

dz 



'dn 



6 qIGM 
°^ " ^ 'dust 



(1 + ^)^ 

a;2 H(z) 



exp 




(27) 



to gain further insight. For small values of x (i.e. 
dust very near the point source), the exponential term 
becomes very large and the curve asymptotically ap- 
proaches zero. The effective scattering angle determines 
where the curve will turn over and thus which region of 
redshift which will contribute most to the integral. Pho- 
tons that are able to scatter through larger angles will 
allow sensitivity that extends to higher redshift gas. 

Figure [5] shows Equation [27] for different combinations 
of a, a, E, and Zg. We choose a — 15", E = 1.5 keV, 
a — 0.5 /im, and Zs = 4 as fiducial values, around which 
each parameter was varied. The curves reach a maxi- 
mum at intermediate redshift, but are often asymmetric. 
The vertical dotted lines demarcate the central region for 
which the area under the curve equals 2/3 of the total 
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Fig. 5. — The curve -r- -tq a-s it changes around the case a. = 15", a = 0.5 /im, E = 1.5 keV, and Zs = 4.0. The vertical dotted lines 



mark the central region that encloses 2/3 of the total area under the curve. 



integrated brightness at the respective observation angle. 
Light that scatters close to the point source must gen- 
erally go through large scattering angles to be observed. 
Therefore, scattered light observed at small angular dis- 
tances from the point source center predominantly traces 
gas closer to the source (Figure [5^). For the same reason 
- low photon energies, which can scatter through larger 
angles, will be able to probe gas closer to the point-source 
(Figure[5lD). The curves in Figure[5]3 are degenerate with 

Figure [5^ because the width of the Gaussian, a oc a/E. 
Large dust grains, which are more strongly forward scat- 
tering than small dust grains, will only contribute ob- 
servable scattering from regions of low redshift. This is 
because gas close to the observer receives light rays that 
are nearly parallel to the line between the observer and 
point-source, making the required scattering angle for 
observation lower compared to dust at intermediate dis- 
tances (Figure [St). Finally, viewing a source at smaller 
redshift will probe gas closer the observer. However, the 
shift in region sensitivity is small; scattering halos around 
objects at 2: = 2 — 4 are mostly sensitive to dust in the 
z = 0.5 — 2 range (Figure [SJi). It is interesting that reduc- 
ing the point source's redshift (z^) does not dramatically 
change the redshift range to which the halo image is sen- 
sitive. A campaign that uses image stacking to search for 



an X-ray scattering halo could therefore use broad bins 
in z. 

6.1. Constraints on IGM homogeneAty 

Figure |4] can also be used to place constraints on the 
homogeneity of the IGM. If a scattering halo is asym- 
metrical on an annulus ^ 1 arcminute in size, the dust 
distribution must be inhomogeneous on ^ 300 kpc scale. 
A symmetric scattering halo would thus indicate a pop- 
ulation of dust grains that is smoothly distributed. Ob- 
servation angles around 10" correspond to spatial homo- 
geneity ^ 100 — 200 kpc in diameter. A scattering halo 
that appears asymmetrical across a 10" annulus would 
thus indicate that the dust being probed is confined to 
the gaseous regions of a galactic halo, potentially part of 
an outflow. 

To illustrate, galactic winds ^ 100 km/s would travel 
about 100 kpc per Gyr. As a first approximation (and 
upper limit), the distance traveled by a dust particle that 
has escaped the potential well of its parent galaxy (at Z2) 
and is thus not subject to any appreciable force of gravity, 
is 

(1 +Th{z) ^''^ 

where v is the terminal velocity of the particle. Assuming 
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V ~ 100 km/s for particles that have been ejected from 
a smaU galaxy at z = 2, a dust grain would travel r^ 200 
kpc by z = 1. 

Note that if a 200 kpc thick shell of ejected dust con- 
tinued to expand freely, it would be ~ 1 Mpc in radius 
by z = and about 75 times less dense than it was was 
at z = 1. A comparable amount of dust grains uniformly 
distributed throughout the IGM would become less dense 
by a factor of 8 between z — 1 and z = 0. This means 
that if a uniformly distributed population of intergalactic 
dust is observed, it is more likely to have been distributed 
during the epoch z > 3 than from z < 2 galaxy outflows, 
which dominate the region R < 100 kpc around them. 

6.2. Constraints on quasar variability 

Scattered photons must travel a slightly longer dis- 
tance than non-scattered photons before reaching the ob- 
server. For small angles, this path difference: 
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requires knowledge about the distance to the dust grain. 
We showed that below the value x — a/2(T, almost no 
light will be scattered into the observer's line of sight. 
Substituting this into Equation [221 and recognizing that 
the maximum path difference occurs for observation an- 
gles ~ cr, we get 
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(30) 
for choices of E and a that give large a. 

The extra distance traveled by a scattered photon leads 
to a time delay At ~ 10 years, which is much shorter 
than t he suggested quasar lifetime, ~ 10^ — 10^ years 
(jMartini 2004) . It is therefore unlikely for an X-ray scat- 
tering halo to be truncated due to a particular quasar's 
activity switching "on" recently. However, dramatic vari- 
ations in a quasar's brightness on < 10'* year timescales 
would become apparent if a scattering halo exhibit radial 
variation. If we traced the region of observable scattering 
from a short burst of X-ray light, it would take the shape 
of an elliptical shell. Thus a ring of increased (or de- 
creased) surface brightness might appear around quasar, 
offering a means to co nstrain the long-term variability of 
AGN activity (e.g. Mi ralda-Escud c 1999). 

7. COMPARISON TO PETRIC ET AL. (2006) 

High spatial angular resolution coupled with a focus- 
ing power that confines ~ 99% of point-source light to a 
2" (4 pixel) radius makes Chandra the best X-ray obser- 
vatory in operation today for observing scattering halos. 
The relatively narrow point-source function (PSF) as- 
sures that there is enough contrast between the source 
and background to observe an accumulation of dust- 
scattered photons a small angular distance away. As 
illustrated below, the most difficult part about observ- 
ing scattering halos due to intergalactic dust is obtaining 
a statistically significant number of photon counts with 
r espect to the niodel s. 

iPetric et al.l (|2006f) used the non-detection of an X- 
ray halo image around QSO 1508-1-5714, a bright quasar 
~ 10^*^ erg/s/cm^ at z = 4.3, to place the limit fljg^ < 
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Fig. 6.— (Top:) The radial profile for QSO 15084-5714 in the 
1-8 keV range, with the model point-spread function in grey. 
(Bottom:) The residual intensity after subtracting the PSF model 
of Gactz ct al. (2004), compared to the intensity expected from a 
p = 4 power-law distribution of grey dust (dashed line) and 1 fim 
size grains (dotted line) with total mass density f^^^^ = 10~^. 

2 X 10~^. They did this by comparing the expected halo 
brightness from a dust population of 1.0 fim size grains 
to the brightness of the PSF wings predicted by MARX 
(a ray-trace simulation) and the wings from a previous 
Chandra observation of 3C 273. The number of scattered 
photons expected from 1.0 fim dust is several times that 
expected from a distribution of dust grain sizes, allowing 
us to relax their constraints. 

The reliability of MARX to simulate the Chandra PSF 
wings is uncertain. High energy X-rays are more subject 
to deviations in scattering by micro-roughness in the high 
resolution mirror assembly (HRMA) and require smaller 
grazing incidence angles that are difficult to obtain while 
testing the HRMA on the ground. In addition, appro- 
priately determining the relative brightness of the PSF 
wings from objects like 3C 273 is difficult due to pile-up 
- that is, the tendency of the Chandra CCD instrumen- 
tation to register an incorrect photon energy or mistake 
multiple photon events for a cosmic ray - which occurs 
for bright objects (> 1 X-ray photon per second). Pile- 
up causes an obvious depression of counts in the center 
of 3C 273, prevents the accurate measurement of source 
flux, and thereby degrades the accuracy of the radial pro- 



TABLE 4 
Observed and expected properties for QSO 1508H 
ITS scattering halo. 



5714 AND 



Energy (keV) 


^obs(< 2"L 


Af=S^(10-500"i^ 


JVbk^ 


0.3- 1.0 


1943 




375 (38)i 


24809 


1.0- 2.0 


1915 




32 


15074 


2.0- 3.0 


637 




3 


12247 


3.0- 4.0 


368 




1 


8479 



^Number of counts observed within a 2" radius of point source 
center. 

''The total number of scattered counts expected to be observed in 
the region 10" — 500" away from the point source, rounded to the 
nearest integer. 

"^The estimated number of background counts, extracted from 
a circular region on the same image, for an area covered by an 
annulus 10 ~ 500" wide. 

"^With a factor of 10 reduction due to absorption. 
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file within - 10", or 20 pixels ()Gaetz et al.ll2004D . Since 
the radial profile from 3C 273 appears to match QSO 
1508+5714 best around 5" (10 pixels), the extrapolated 
wing profile is more likely to be an overestimate of the 
true PSF - perhaps explaining why the intensity from 3C 
273 matches the upper limit s on the instrume ntal mea- 
surement between 30 - 60" (jPetric et al .l l2006l Fig. 21. 

We use the PSF calibration model by iGaetz et al.l 
()2004[ ) to gauge the PSF + halo profile that would have 
resulted from the p = 4 dust distribution with density 
n]GM ^ iQ-5 A quick examination of QSO 1508+5714 
(obsid 2241) using a circular region with a 4 pixel ra- 
dius yields approximately 5400 counts in the 0.3-8 keV 
range. (This is a lower limit to the total number of source 
counts, as it does not consider counts spread into the 
PSF.) For several energy bins. Tabic 0] shows the number 
of point source counts, expected number of halo counts 
in an annulus 10 — 500" centered on the point source, and 
the number of background counts expected in that same 
region. The majority of the scattered light will have pho- 
ton energies < 2 keV. However, for soft X-rays (< 1 keV) 
and large dust grains, the RG cross-section can overes- 
timate the intensity of scattering by a factor ^ 10 du e 
to absorption by grain material ([Smith fc Dwekl [19981 ). 
For this reason, the 375 scattered counts expected in the 
10 — 500" annular region are likely too large by an or- 
der of magnitude. We therefore exclude photon energies 
below 1 keV from further evaluation. To take advantage 
of the soft X-ray band, the halo profile must be modeled 
with Mie scattering. 

Figure ini compares the radial brightness profile of QSO 
1508+5714 in the 1-8 keV energy band to the PSF and 
dust models. The mag nitude of th e erro r bars show 
that the observations of IPetric et al.l (|2006f) are consis- 
tent with iljj^^ < 10~^ for the power-law distribution 
of dust grains chosen. However, the brightness profile is 
also consistent with a population of IGM dust composed 
purely of 1 /im sized grains. By inclu ding energies below 
1 keV and using the RG cross-section. IPetric et al.l ()2006f ) 
were expecting a much brighter X-ray halo. This was an 
overestimate leading to the limit ^]^^ < 2 x 10^®, which 
our analysis shows can be relaxed to ^^^^ < 10^^. This 
emphasizes the need for a more detailed modeling, incor- 
porating statistics and Mie scattering, to be implemented 
in the future. An even more important task would be the 
careful calibration of the Chandra HRMA PSF between 
1 and 20". 

8. DISCUSSION 

X-ray scattering provides a unique view of the cosmic 
history of star formation, feedback, and the IGM that is 
complementary to UV, optical, and infrared studies. A 
summary of our main results are as follows: 

(i) If dust is evenly distributed throughout the IGM, 
with n]^^ = 10"*^, the dust-scattered X-ray light will be 

~ 5% (-B^cv) ^^ ^^^ point source brightness for objects 
at z > 2. The intensity of the scattered light is directly 
proportional to fijju^. 

(ii) Also in the uniform IGM case, point sources at 
larger values of z trade larger dust column densities for 
a narrower halo profile, making it slightly more difficult 
to distinguish scattered light from the point source. 

(iii) For high-z point sources, scattered light comes 



TABLE 5 

Systematic offset in magnitude for an intergalactic 

population of grey graphitic grains. 





(P = 


= 3) 


(P = 


= 4) 


z 


Av 


-4i.78 


Av 


Ai.78 


0.50 


0.008 


0.007 


0.012 


0.010 


1.00 


0.010 


0.010 


0.017 


0.014 


2.00 


0.012 


0.012 


0.019 


0.017 


4.00 


0.013 


0.013 


0.021 


0.019 



predominantly from intergalactic dust at z < 2. If the 
majority of IGM dust came from star formation within 
sub-L* galaxies at z > 3, assuming a constant comoving 
density of dust grains may still be reasonable for point 
sources at z = 4. 

(iv) If X-ray light scatters through a dense, dusty 
clump in intergalactic space, the resulting halo image 
will have a fiatter radial profile compared to the uni- 
formly distributed case. The further away the clump is 
from the background point-source, the larger the halo 
image will appear to an observer. 

(v) If the dust mass to metallicity ratio is relatively 
constant, as observed locally, then a typical DLA will 
have an optical depth to 1 keV X-ray scattering ^ 1 — 3%, 
proportional to (1 + Zg)~^. This result scales propor- 
tionally with metallicity, so DLAs with near-solar abun- 
dances may have Tx > 5%. 

(vi) Dusty clumps at z > 1 would have to be > 200 kpc 
in size for a halo image to appear at angles larger than 
10" from the point source center. X-ray scattering thus 
offers the opportunity to test if dusty outflows exist near 
quasar absorption systems, but the outflow may require 
column densities of dust that are so far unprecedented. 
Perhaps an X-ray point source whose line of sight is close 
to a foreground galaxy can test the ability of large grains 
to be expelled from L* galaxies. 

(vii) A scattering halo that is symmetric on an annulus 
of radius > 60", would indicate that IGM dust is dis- 
tributed uniformly. A scattering halo that is asymmetric 
on an annulus of radius 10" would indicate the the dust 
distribution is traced mainly by galaxies {R < 50 kpc) 
and is not distributed very far out in the galactic halo 
potential or the IGM. Either result would provide clues 
to galactic and star formation feedback mechanisms. 

(viii) The previous constraint of ^^^^ < 2 x 10^^ must 
be relaxed when considering the absorption of X-ray light 
below 1 keV. Photons observed at 0.5 keV today will be 
shifted above 1 keV for z > 1, but in Section[n]we showed 
that the scattering halo is more sensitive to material at 
z < 2 and large dust grains at low redshift. Therefore, 
what would be the dominant contribution to the scatter- 
ing calculation, X-rays < 1 keV and 1 /im sized grains, is 
overestimated by Rayleigh-Gans scattering. Deeper ob- 
servation of high redshift quasars, in conjunction with a 
more careful calculation involving the Mie scattering so- 
lution, is required to place a stronger constraint on the 
intergalactic dust distribution. 

To conclude, we briefly note that the dust models used 
throughout this paper would create meaningful system- 
atic errors for both optical and infrared extragalactic sur- 
veys. Intergalactic dust dims the most distant objects, 
and if not accounted for, would affect cosmologists' abil- 
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ity to place more stringent constraints on dark energy 
models. A systematic shift as low as Sm ~ 0.02 would 
reduce the effectiveness of future supernovae surveys to 
constrain the dark energy equation of state parameter w 
lyircv & Ealct 2007). A systematic shift in magnitude 
roughly co rrelates as Sm ^ —0.5 dw and dm ^ — SQm for 
z < 1 fe.g. iZhandbOOa iMenard et al.ll2QT0al: iCorasanitil 
|2006[) . Thus an extinction offset ^ 1% would create an 
offset comparable to the error bars sought after in preci- 
sion cosmology. 

A distribution of graphitic grains larger than 0.1 /im is 
grey even to the reddest filters (0.893 /xm and 1.78 /im) 
to be used in cosmological supernovae surveys (LSST 
and WFIRST, respectively). We calculate the systematic 
offset in magnitude for our grey graphitic grain models 
using n^^^ = 10^5 (Table [5]). There is little offset be- 
tween Ay and A1.78 - the extinction at 1.78 /im - which 
is th e center of the red dest filter on the WFIRST mis- 
sion (jGreen et al.ll201ll ) . For supernovae surveys beyond 
z > 1, 5m, ~ 0.01 — 0.02 in both the visual and infrared 
bands. Such small systematic offsets would call for a 
deep field survey, rather than wide field survey that in- 
creases the precision of t he photometric mea surements 
but not their accuracy (|Virev fc Ealetl l2007t ) . Obser- 



vational searches for X-ray scattering halos due to in- 
tergalactic dust would constrain these systematic offsets 
and help determine the best future course of action for 
large-scale optical and infrared cosmological surveys. 

Finally, we note that the values in Table [5] differ by 
those calculated by Aguirre (; 19991 ) for a similar cosmic 
density of intergalactic dust. For the model p — A 
graphite grain distribution, kv = 2.4 x 10^ cm^/g. 
Using this val ue with ^^^ = 10^^ in Equation 9 of 
lAguirrd (|1999[ ) implies Av{z = 0.5) - 0.02. However, 
Aguirre uses Vim = 0.2 and Ha — Q- In that case, the 
distance (and hence the column density) one would 
calculate for z = 0.5 is larger by a factor of two than 
that in the currently accepted ACDM universe. This 
brings previous work into agreement with the results 
shown in Table [5] 

We thank the anonymous referee for their attention to 
detail and help in greatly improving the clarity of the 
paper. We thank J. G. Peek and J. Schroeder for useful 
discussion about extragalactic dust and cosmology. This 
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der the NASA Earth and Space Science Fellowship Pro- 
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